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Abstract 

The method of self-similar factor approximants is applied to calculating the critical 
exponents of the 0(A^)-symmetric (/?^ theory and of the Ising glass. It is demonstrated 
that this method, being much simpler than other known techniques of series summa- 
tion in calculating the critical exponents, at the same time, yields the results that 
are in very good agreement with those of other rather complicated numerical meth- 
ods. The principal advantage of the method of self-similar factor approximants is the 
combination of its extraordinary simplicity and high accuracy. 



PACS: 

05.70.Jk Critical point phenomena in thermodynamics, 
02.30.Lt Sequences, series, and summability, 
02.30.Mv Approximations and expansions 

1 Introduction 



The knowledge of critical exponents, characterizing critical phenomena, provides us with ba- 
sic information on the behavior of thermodynamic quantities in the vicinity of critical points 
[1-3]. This is why so much efforts have been devoted to the experimental measurements as 
well as to theoretical calculations of these exponents. 

Because of the complexity of realistic theoretical models, critical exponents can usually 
be represented by power series obtained with the help of some perturbation theory. Such se- 
ries are practically always divergent, which requires to use resummation techniques allowing 
for the determination of effective limits of divergent series. The standard approach applied 
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to the summation of series, associated with critical indices, is based on the Pade-Borel-Leroy 
transformation and its variants [4]. Another approach is based on optimized perturbation 
theory [5], where the resummation is due to control functions transforming divergent series 
into converging ones. Introducing control functions through the change of expansion vari- 
ables, it is possible to resum the field-theoretic expansions for critical exponents [6,7]. The 
convergence of the optimized perturbation theory can be greatly accelerated by invoking the 
self-similar approximation theory [8-16], as has been done for calculating critical exponents 
[17]. All these approaches, mentioned above, require quite complicated numerical calcula- 
tions. A purely numerical procedure of calculating critical exponents is due to Monte Carlo 
simulations [18-29]. 

In the present paper, we suggest a new approach for the summation of series related 
to critical exponents. This approach uses the method of self-similar factor approximants 
[30-33], whose mathematical foundation is based on the self-similar approximation theory 
[8-16]. The construction of the factor approximants is very simple and straightforward. We 
recall the main definitions in Section 2. Then, in Section 3, we apply these approximants 
for the summation of the ^-expansions for the critical exponents of the A^-vector (p^ field 
theory. Our very simple method yields the results that arc in perfect agreement with the 
most complicated numerical procedures. In Section 4, we demonstrate that the suggested 
method is applicable even for such a notoriously difficult problem as finding the critical 
indices from the high-temperature series expansions for spin glasses. Finally, Section 5 is 
conclusion. 



2 Self-similar factor approximants 

Suppose that our aim is to reconstruct a real function f{x) of a real variable x, when the 
function is represented by its asymptotic expansion at a; as a power series 

k 

fk{x)^J2anx'', (1) 

n=0 

where A; = 0, 1, 2, . . .. Without the loss of generality, we may assume that oq = 1 in expansion 
(1). This is because if instead of form (1) we would have a more general expression 

n=0 

with a given function then we could return to expansion (1), with ao = 1, by defining 

The sequence {fki^)} is usually divergent for any finite x. 

The method of self-similar factor approximants [30-33] makes it possible to extrapolate 
the asymptotic expansion (1), valid only for a; — > 0, to the whole region of the variable a; > 0. 
When k = 2p is even, with p being an integer, then the even-order factor approximant is 

fUx)^f{{l + A,xr . (2) 
1=1 
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The parameters Ai and Ui are obtained from the re-expansion procedure, that is, by expand- 
ing approximant (2) in powers of x up to the k-th order and equating the latter expansion 
with the initial one given by form (1). This re-expansion procedure yields the set of 2p 
equations 

Y,niA^ = Bn (n=l,2,...,2p) (3) 

i=l 

for 2p parameters rii and Ai, with the right-hand sides 

Bn = -TT lim - — In fhix) . 

As is evident, the quantities n^, Ai, as well as depend on the considered order k. But for 
avoiding excessively cumbersome notations, we do not use here the double labelling. When 
/c = 2p + 1 is odd, the odd-term factor approximant is 

/2p-M(^) = n(i+^^^r , (4) 

with the parameters Ai and Ui defined by the re-expansion procedure yielding the set of 
equations 

p+i 

Y^mA^i^B^ (n = l,2,...,2p+l), (5) 

with the scaling condition Ax = \. 

In this way, for any given expansion (1), the construction of self-similar factor approx- 
imants (2) or (4) is rather simple and straightforward. It has been shown [30-33] that the 
factor approximants are more general and accurate than Pade approximants, having, in ad- 
dition, a principal advantage of being uniquely defined. This means that for each given order 
k of expansion (1) there is just the sole factor approximant, while for each k there exists a 
table of k different Pade approximants P[M/iv], with M + N = k. There is no general recipe 
for choosing one of the k available Pade approximants. One often chooses the diagonal ones, 
but, as is easy to show, the latter are not always the most accurate ones. Such a problem 
of multiple possibilities does not arise for factor approximants: for each k, of the expansion 
fk{x), there is just one factor approximant fk{x). 



3 Exponents for 0(A^)-symmetric theory 

Let us consider the A'"-component vector </7^ field theory, for which the critical exponents 
can be obtained in the form of the ^-expansions, with e = 4 — d, and d being the space 
dimensionality. The derivation of these dimensional expansions can be found in the book 
[6]. In the Appendix A. wc give the expansions that are considered in the present section. As 
is known, such expansions are divergent and require a resummation procedure. To illustrate 
more explicitly how the method of self-similar factor approximants works, let us start with 
the scalar single-component field (N — 1). Then we have 

T] ~ O.OlSSe^ + 0.0187£^ - 0.0083^^ + 0.0257^^ , 
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i/-^ ~ 2 - 0.333£ - O.lUe^ + O.UAs^ - O.SOTs^ + 0.95l£^ , 

ujc^e- 0.63£^ + 1.62£^ - 5.24£^ + 20.75e^ . (6) 
We reduce each of these expansions to the form 

k 

^(£)=/o(£)^a„£^ (7) 

n=0 

where Qq — 1. According to Section 2, we construct the factor approximants 

fkie)-^Me)Ui^ + ^i^r , (8) 

i=l 



in which 



|, A: = 2p = 2,4,... 
%i , A; = 2p + 1 = 3, 5, . . . 



and the parameters Ai and rii arc obtained from the re-expansion procedure. Setting s — 1, 
we find the desired approximation /^(l) = f^- The error bar for the approximant is given 
by 

±^(/:-/:-i) ^ = 2,3,... 

Reducing the series for r] to form (7), we have 770 (s) = 0.0185£^. Constructing the factor 
approximants (8), we find for (^) the parameters 

Ai = 1.898511 , ni = 0.532423 , 

and for ri^{e), we get 

Ai^l, ni = 0.789878 , A2 = 5.110862 , 713 = 0.043228 . 
Setting £ = 1, we obtain 

7]; = 0.032602 , ri; = 0.034588 . 
Thus, we conclude that the factor approximants give 

77 = 0.035 ± 0.001 . 

In the same way, we proceed with the series for In the second order, we find 

Ai = -0.869203 , = -0.191555 . 

For the third order, we have 

Ai = 1 , Til = -0.152346 , A2 = 0.023910 , 713 = 13.335389 . 
In the fourth order, we get 

= 3.027805 , 7ii = -0.006791 , A2 = -0.440821 , 712 = -0.424352 . 
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And for the fifth order, we find 

Ai = l, m = -0.045336 , A2 = 4.168053 , = -0.001772 , 

A3 = -0.312951 , rig = -0.700494 . 
Setting £ = 1, for the factor approximants ul, we obtain 

ly* = 0.738227 , i/* = 0.616528 , = 0.633852 , i/* = 0.628417 . 

Hence, the result is 

ly = 0.628 ± 0.003 . 

Finally, reducing the series for cu to form (7), we have u!o{s) — £■ Following the standard 
procedure, for the factor approximant uj2{s), we get 

Ai = 4.512857 , rii = -0.139601 . 

For ou^{e), we find 

Ai^l, m = 0.006238 , A2 = 4.547986 , ris = -0.137151 . 
And for ujI (e) , it follows 

Ai = 4.511659 , ni = -0.139637 , A2 = 107.494872 , 712 = -0.7 x 10"^ . 
Setting £ = 1, we obtain 

uj^ = 0.787958 , uj^ = 0.787160 , ujI = 0.787934 . 
Hence, we come to the value 

u = 0.788 ±0.0004 . 
Other critical exponents can be obtained from the scaling relations 

a^2-ud, /3=^(rf-2 + r/) , 7 = z/(2 - r/) , 5 = ^±|^ , (9) 

which for the dimensionality d — 3 simplifies to 

a^2-3u, /?=^(l + ry), ^^^2-7?), <^ = ^ ■ (10) 

2 1 + 77 

Using here the found results for the factor approximants, we have 

a = 0.116 ±0.009, /3 = 0.325 ± 0.002 , 7 = 1.234 ± 0.005 , 5 = 4.797 ± 0.006 . 

Accomplishing in the same way calculations for the arbitrary number of components N, 
we obtain the factor approximants for the critical exponents using the general series from 
the Appendix A. Our results are presented in Table 1. It is worth emphasizing that in the 
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two limiting cases of iV = —2 and N = oo, where the exact critical exponents are known, 
our results coincide with these exact values. For N — —2, the exact exponents are 

in any dimension. And in the limit of large N, the exact exponents are 

d-A ^1 2 , d+A 

77 = 0, T^=-j^, uj^A-d (N^oo), 
where d is dimensionality. In three dimensions, the latter transforms to 

« = /3=^, 7 = 2, 5 = 5, r? = 0, a; = l. 

Since our results tend to the exact values when N oo, the error bars diminish for ^ 1, 
tending to zero, as iV — > oo. Thus, for N = 100, the error bar is 10~^, for N = 1000 it is 
10~^, and for N — 4:, the error bar is 10~^. The error bars for A" 10 diminish as 1/A". 
The error bars for the factor approximants, up to A" = 10, obtained from the expansions for 
T), u, and uj, are shown in Table 2. 

Critical exponents have been calculated by Monte Carlo simulations [18-29,34-39] and 
other complicated numerical methods, as is reviewed in Ref. [6,17,40-42]. Our results in 
Table 1 are in very good agreement with all these calculations. The advantage of our method 
is its simplicity. We have used only the expansions from the Appendix A. We do not need 
to know the large-order behavior of £-expansions, which is required for other methods. 



4 Exponents for spin glass 

Here we show that the method of self-similar factor approximants can be applied to such a 
notoriously difficult problem as summing the high-temperature series of the Ising spin glass. 
This model is described by the Hamiltonian 

// = - ^ JyCTifJj- , (11) 

(ij) 

in which (ij) implies the summation over nearest neighbors, Ci takes values ±1, and Jij are 
independent random variables, whose dimensionless forms J^- = pJij, with P being inverse 
temperature, occur with the probability 

p(Jij) = I [SiJij -J) + S(Jij + J)] (12) 

where J^j = f3J is a parameter. Monte Carlo simulations [43-45] demonstrate the existence 
of the phase transition in three dimensions. 
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The phase transition corresponds to a singularity in susceptibihties. One considers two 
types of the latter, the Edwards- Anderson susceptibility 

Xea = ^ E < (< ^i^3 >') > ' (13) 
and the auxiliary susceptibility 

i,3 

Here is the total number of lattice sites, single angular brackets < . . . > refer to thermal 
averaging, and the double brackets <^ . . . ^ refer to averaging with respect to the distri- 
bution of interactions, defined by probability (12). When temperature T approaches the 
critical temperature Tc, susceptibilities (13) and (14) behave as 

XEA « (T - Te)-^ , x' oc (T - %)-^' . (15) 

High-temperature series expansions for susceptibilities (13) and (14) are represented as 
series in powers of 

w = tanh2(/3J) . (16) 

Analyzing the series 

XEA ^ E ^nW"" , X' - E ' (1 '') 

n n 

one aims at finding the critical exponents 7 and 7' characterizing the critical behavior (15). 
These exponents are connected with each other through the scahng relations 

7 = (2-r/)i/, 7' = (4-d-2r/)i/, (18) 

where 1] and u are the critical exponent defining the behavior of the correlation function 
<C (< aiaj >^) ^ and the correlation length ^ oc (T — Tg)^'', and where d is dimensionality. 
Therefore, the exponents 7 and 7' can be expressed one through another by means of the 
scaling relation 

27 = 7' + i/ci, (19) 

provided i/ is known. 

The analysis of expansions (17) turned out to be extremely difficult. This is because 
the first few terms of the series contain little information on spin-glass ordering. Actually, 
the coefficients of the first three terms of the series for xea are identical to those for the 
susceptibility series of the pure Ising model. In fact, one cannot see any spin- glass behavior 
until one gets contribution from higher orders. This happens because an essential feature of 
spin glass is frustration, which reveals itself only in higher orders of the series. Hence, any 
analysis, depending sensitively on the first few terms in determining the critical exponents 7 
and 7', is not hkely to give correct answers. These difficulties have been described in detail 
by Singh and Chakravarty [46,47], who found that more than ten terms in expansions (13) 
and (14) are necessary to be able to estimate the critical exponents 7 and 7'. They derived 
[46] in three dimensions expansions for Xea and x' up to 17-th order and in four dimensions, 
an expansion for xea up to 15-th order. 
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However, even having quite a number of terms in expansions (17), it is very difficult to 
find the related critical exponents. Since, as is stressed above, the effects of frustration reveal 
themselves only in high orders of expansions, so that the lower orders do not provide correct 
information on spin glass behavior. The method that has been found [47] to be most suitable 
to this problem is that of inhomogeneous differential approximants, which is a generalization 
of the d-\og Pade summation. A weak point of this method is that the approximants, for 
each given order k of an expansion, are not uniquely defined. Thus, for an expansion of order 
k = 10, there are 42 variants of the approximants, for the expansion of order k = 15, there 
are 96 variants, and for the order k = 17, one has 136 variants of different approximants. 
All these various approximants yield the results that are quite different from each other, 
and it is not clear which of them are to be accepted as correct and which as wrong, so that 
the problem arises of a subjective, not strictly defined, choice of some of them labelled as 
"well-behaved" . 

Now let us apply to expansions (17) the method of self-similar factor approximants of 
Section 2. In three dimensions {d = 3), both series for xea as well as for x' are known [47] 
up to seventeenth order. The effect of frustration, typical of spin glass, occurs in the series 
for x' much earlier than in that for xea- Therefore more accurate results can be obtained 
considering the series for which are 

X' = 1 + + W2w^ - 192w^ + mSw^ - 758W + 42822^^ - 221856^^^+ 

+ 1147878w^° - 5980608w;^^ + 32318910?^/^ - 167464128^^^ + 906131742^^^- 

-4849958304^^^ + 25952889798w^^ - 141648771168w;^^ . (20) 

For these series, we construct the factor approximants following the standard procedure of 
Section 2. The closest singularity to the origin defines the critical points Wc and the related 
critical exponents 7'. For the series of order k — 15, we find Wc — 0.42 and 7' = 2.07; for 
k = 16, we have Wc = 0.39 and 7' = 1.44; and for k = 17, we get Wc = 0.41 with 7' = 1.82. 
Thus, for the critical exponent 7', we obtain 

7' =1.82 ±0.19 {d = 3). (21) 

This can be compared with the results of the Monte Carlo simulations [45], which find the 
phase transition at = 1.18 ± 0.03, with the critical exponent 

7;,c; = 1-87 ±0.28 (d = 3). (22) 

Using our result (21), the known value [40] oi u — 1.3 ± 0.1, and the scahng relation (19), 
we find the critical exponent 

7 = 2.86 ± 0.24 (d = 3) . (23) 

From the Monte Carlo simulations [40] it follows 

7mc = 2.89 ± 0.29 (d = 3) . (24) 

As is seen, the critical exponents in Eqs. (21) and (23) are close to the Monte Carlo values 
in Eqs. (32) and (24), respectively. 
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In four dimensions (d = 4), only the series of fifteenth order for xea are available [47], 
which are 

X = 1 + 8w + 56w^ + 392w^ + 2408w^ + 15272w^ + 85352w^ + 508808^ V 

+2625896w^ + 15111976w^ + 72067672w^° + 421464680w^^+ 

+1851603192^^^ + 11810583208^^^ + 46346625320^^^ + 347729503368^/;^^ . (25) 

Constructing the factor approximants for these series, in the highest orders we find = 0.20, 
with 7 = 1.59, for k = lA and Wc — 0.21, with 7 = 2.35, for k — 15. Therefore, for the 
critical exponent 7, we obtain 

7 = 2.35 ± 0.38 (c? = 4) . 

To our knowledge, Monte Carlo simulations for d = 4 are not available. And the method of 
inhomogeneous differential approximants [47] estimates 7 ~ 2.0 ±0.4. Since there is neither 
an expansion for x' iior information on other indices, it is not possible to determine the 
exponent 7' in 0? = 4. 

The example of the present Section shows that the method of self-similar factor approx- 
imants can be applied to rather complicated series with very nontrivial behavior, requiring 
the consideration of high-order terms. 

5 Conclusion 

The method of self-similar factor approximants [30-33] is apphed to calculating the critical 
exponents of the A^^-component vector (^^ field theory and of the Ising spin glass. The first 
example is chosen because of the wide interest to the 0(A^)-symmetric ip"^ theory, which 
serves as a typical model for characterizing the critical behavior of a large variety of physical 
systems. We showed that the apphcation of the method to £-expansions is very simple 
and straightforward at the same time providing the accuracy comparable with that of other 
essentially more complicated techniques. 

The case of high-temperature expansions for the Ising spin glass is taken as an example 
of series with a notoriously nontrivial structure, requiring the consideration of high-order 
terms and making it very difficult an unambiguous determination of the exponents for sus- 
ceptibilities by other known methods. Our method allows us to find the exponents that are 
in good agreement with Monte Carlo simulations, when the latter arc available. 

In the present paper, we have concentrated on the calculations of critical exponents. Of 
course, determining the critical points is also of importance. For instance, recently there has 
been a great interest to an accurate calculation of the critical temperature Tc for interacting 
Bose gas (see review articles [48,49]). The most accurate results have been obtained so 
far by using the ideas of the optimized perturbation theory [5] in Refs. [50-55] and by 
employing Monte Carlo simulations [56-60]. These are rather involved numerical techniques. 
The method of self-similar factor approximants can also be applied to this problem, which, 
however, is a topic for a separate investigation. 
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Appendix A 



The derivation of the expansions for the A^-component field theory, considered in Section 
2, can be found in the book [6]. These expansions are 



{N + 2)e^ 

2(iV + 8)2 ' 4(Ar + 8)2 



1 + 



-A^2 + 56A^ + 272]- 



16(iV + 8)4 

^3 



bN^ + 230A^^ - 1124Ar2 _ 17920A^ - 46144 + 384C(3)(iV + 8)(5iV + 22) 



ISA^e + 946Ar5 + 27620A^^ + 121472A^^ - 262528A^^ - 2912768A^ - 5655552- 



64(Ar + 8)6 

-16C(3)(A^ + 8) (n^ + lOA^^ + 1220A^^ - lUm'^ - 68672A^ - 171264) + 

+1152C(4)(iV + 8)^(5iV + 22) - 5120C(5)(iV + Sf{2N'^ + 55A^ + 186)] } , 



i/-^ = 2 + 



{N + 2)e 
N + 8 



2{N + 8Y 



[13N + 44]+ 



r-3 



3N^ - 452N^ - 2Q72N - 5312 + 96C(3)(iV + 8)(5iV + 22) 



+ 



+ 



3A^^ + 398A^^ - 12900A^'^ - 81552A^2 - 219968A^ - 357120+ 



8(Ar + 8)6 

+16C(3)(A^ + 8) (3A^^ - 194Ar3 + u^n"^ + 9472A^ + 19488) + 288C(4)(A^ + 8)^(5A^ + 22)- 



-1280C(5)(iV + 8f (2Ar2 + 55Ar + 186) 



+ 



3A^^ - 1198A^^ - 27484A^^ - 1055344A^^ - 52421 12 A^^ - 5256704A^V 



128(iV + 8)8 

+69990407V-626688-16C(3)(7V+8) (l37V^ - 310Ar^ + 19004Ar^ + 1024007V2 - 3815367V2- 
-2792576A^ - 4240640) - 1024(^(3) (A^ + 8)^ (2A^^ + 18A^^ + 981 A^^ + 6994A^ + 11688) + 
48C(4)(A^ + 8)-'^ (SA^^ - 194A^^ + lA&N"^ + 9472iV + 19488) + 
+256C(5)(Ar + 8f {ibbN^ + 3026Ar3 + 989Ar2 _ 56018A^ - 130608) - 
-6400C(6)(Ar + 8)^ (2A^^ + 55A^ + 186) + 56448C(7)(Ar + 8)^ (mA^^ + 189A^ + 256)] } , 



uj{e) — e — 



3e' 



{N + 8y 



[3Ar + 14] + 



"4(Ar + 8)4 



33 A^^ + 538A^2 + 4288A^ + 9568 + 96C(3)(Ar + 8)(5A^ + 22) 



+ 
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<r4 

5N^ - 1488A^^ - 46616A^^ - 419528A^^- 



16(A^ + 8)6 

-1750080A^ - 2599552 - 96C(3)(A^ + 8) (63A^^ + 548A^^ + 1916A^ + 3872) + 



+288C(4)(A^ + 8)3(5iV + 22) - 1920C(5)(iV + 8)^ (2N^ + 55iV + 186 

13A^^ + 7196A^^ + 240328A^^ + 3760776A^^+ 



+ 



64{N + 8Y 

+38877056Ar3 ^ 223778048Ar2 + 660389888Ar + 752420864- 
-16C(3)(A^+8) (9A^^ - 1104A^5 - 11648A^^ - 243864A^=^ - 2413248A^2 - 9603328A^ - 14734080) 

-768C^(3)(iV + 8)2 (6A^^ + 107A^=^ + 1826A^2 + 9008iV + 8736) - 

-288C(4)(iV + 8)^ (63iV3 + 548Ar2 + 1916Ar + 3872) + 

+256C(5)(A^ + Sf (305N^ + 73S6N^ + 45654Ar2 + 143212A^ + 226992) - 

-9600C(6)(iV + 8)^ (2N^ + 55iV + 186) + 112896C(7)(iV + 8)^ (lAN^ + 189iV + 256) 
Here e — 4 — d is assumed to be asymptotically small, £ ^ 0. 
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Table Captions 



Table 1. Critical exponents for the A'"-component </7^ field theory, obtained by the 
summation of ^-expansions using the method of self-similar factor approximants. 

Table 2. Error bars for the critical exponent of table 1. For large N ^ 10, the error 
bars diminish with N as 1/A^. 
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Table 1 



N 


a 


P 


7 


S 


V 


u 


cu 


-2 


0.5 


0.25 


1 




5 







0.5 


0.80118 


-1 


0.36844 


0.27721 


1.07713 


4 


88558 


0.019441 





54385 


0.79246 





0.24005 


0.30204 


1.15587 


4 


82691 


0.029706 





58665 


0.78832 


1 


0.11465 


0.32509 


1.23517 


4 


79947 


0.034578 





62854 


0.78799 


2 


-0.00625 


0.34653 


1.31320 


4 


78962 


0.036337 





66875 


0.78924 


3 


-0.12063 


0.36629 


1.38805 


4 


78953 


0.036353 





70688 


0.79103 


4 


-0.22663 


0.38425 


1.45813 


4 


79470 


0.035430 





74221 


0.79296 


5 


-0.32290 


0.40033 


1.52230 


4 


80254 


0.034030 





77430 


0.79492 


6 


-0.40877 


0.41448 


1.57982 


4 


81160 


0.032418 





80292 


0.79694 


7 


-0.48420 


0.42676 


1.63068 


4 


82107 


0.030739 





82807 


0.79918 


8 


-0.54969 


0.43730 


1.67508 


4 


83049 


0.029074 





84990 


0.80184 


9 


-0.60606 


0.44627 


1.71352 


4 


83962 


0.027463 





86869 


0.80515 


10 


-0.65432 


0.45386 


1.74661 


4 


84836 


0.025928 





88477 


0.80927 


50 


-0.98766 


0.50182 


1.98402 


4 


95364 


0.007786 





99589 


0.93176 


100 


-0.89650 


0.48334 


1.92981 


4 


99264 


0.001229 





96550 


0.97201 


1000 


-0.99843 


0.49933 


1.99662 


4 


99859 


0.000235 





99843 


0.99807 


10000 


-0.99986 


0.49993 


1.99966 


4 


99986 


0.000024 





99984 


0.99979 


oo 


-1 


0.5 


2 




5 







1 


1 
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Table 2 





T] error bar 


v error bar 


cu error bar 


-2 








0.0280 


-1 


0.0007 


0.0008 


0.0013 





0.0010 


0.0018 


0.0110 


1 


0.0010 


0.0027 


0.0043 


2 


0.0009 


0.0034 


0.0016 


3 


0.0008 


0.0038 


0.0020 


4 


0.0007 


0.0039 


0.0016 


5 


0.0006 


0.0038 


0.0007 


6 


0.0005 


0.0036 


0.0005 


7 


0.0004 


0.0034 


0.0019 


8 


0.0003 


0.0032 


0.0032 


9 


0.0002 


0.0030 


0.0042 


10 


0.0001 


0.0029 


0.0048 
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